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We propose and study an "intrinsic probing" approach, without introducing any external detector, 
to mimic cavity QED effects in a qubit-nanomechanical resonator system. This metallic nanome- 
chanical resonator can act as an intrinsic detector when a weak driving current passes through it. 
The nanomechanical resonator acts as both the cavity and the detector. A cavity QED-like effect 
is demonstrated by the correlation spectrum of the electromotive force between the two ends of 
the nanomechanical resonator. Using the quantum regression theorem and perturbation theory, we 
analytically calculate the correlation spectrum. In the weak driving limit, we study the effect on the 
vacuum Rabi splitting of both the strength of the driving as well as the frequency-detuning between 
the charge qubit and the nanomechanical resonator. Numerical calculations confirm the validity of 
our intrinsic probing approach. 

PACS numbers: 85.85.+j, 85.25.Cp 



I. INTRODUCTION 

Recently, nanomechanical resonators (NAMRs) are 
attracting considerable attention (see, e.g., Refs. [H, 
[2,3,4). Also, mechanical analogues of cavity QED 
have been theoretically studied in coupled systems be- 
tween nanomechanical resonators and superconduct- 
ing qubits (see, e.g., Refs. 5,6,7.8,9). Various effects 
in these nanomechanical QED systems were investi- 
gated, including: quantum measurementsS,, the quan- 
tum squeezing of the NAMRs^ i^°'^^ , and the cooling of 
the NAMRsi^ii^iMii^ii^. Some of these theoretical pro- 
posals have recently become experimentally testable due 
to the recent advances in NAMRs and superconducting 
qubits. Numerous Josephson-junction-based supercon- 
ducting qubits have been experimentally realized (see, 
e.g., the review o^*''-^^!"^^'^" ), while studies on NAMRs with 
vibration frequencies of the order of a GHz are approach- 
ing the quantum regime. 

References2ii2^ recently studied a NAMR coupled to 
a double-quantum dot. In Ref. [2l|, the spectrum of the 
transport current was used to study the quantum be- 
havior of this system. The electron transport through a 
mobile island (i.e., a nanomechanical oscillator) with two 
energy levels was studied in Ref.H^ There, the qubit was 
embedded in the NAMR. 

To study cavity QED analogues in a NAMR-qubit sys- 
tem, a crucial issue is how to make the quantum mea- 
surement on this coupled system. Quantum measure- 
ments involve subtle interactions between the system and 
the detector. To carry out a quantum measurement, an 
external probing instrument is typically coupled to the 
measured system. Examples of this include: a single 
electron transistor— coupled to a charge qubit, a trans- 
mission line resonato r^^i^^i^^ coupled to a charge qubit. 



or a shunted dc-SQUID coupled to a flux qubit (see, e.g., 
Refs. HIHIM). In general, these coupled systems can 
be modeled by the Jaynes-Cummings Hamiltonian and 
demonstrate several analogues to cavity QED effects^i^, 
such as vacuum Rabi splitting. These effects can be used 
to verify the coherent coupling between a superconduct- 
ing qubit and a measuring device. 

This study is mainly motivated by recent experiments 
on a high-frequency metallic NAMI^. Previously, non- 
metallic NAMRs were often studied and therefore no 
efficient current would pass through these non-metallic 
NAMRs, and thus no mechanical force acting on the 
NAMRs could be induced to implement quantum mea- 
surements. In this case, an external instrument needs to 
be integrated to probe the coupling between the NAMR 
and the qubit. Here, we study how to probe a cavity QED 
analogue for a metallic NAMR'^^ coupled to a supercon- 
ducting qubit without introducing an external detector. 

In this proposal, the information on the coherent cou- 
pling between the superconducting charge qubit and the 
metallic NAMR can be read out by measuring the in- 
duced electromotive force between the two ends of the 
NAMR. This electromotive force is generated by a cur- 
rent passing through a metallic NAMR in which a mag- 
netic field is applied. There are at least two advantages 
for this intrinsic probing mechanism: (i) the coupling be- 
tween the metallic NAMR and the qubit can be turned 
on or off by the externally-applied voltage, and then the 
information can be read out in a controllable way; (ii) no 
external probing instrument needs to be introduced, in 
contrast to the proposal in Ref. Is^ . 

This paper is organized as follows. In Sec. II, we de- 
scribe the proposed model, and write the Hamiltonian for 
a charge qubit interacting with a driven metallic NAMR. 
In Sec. HI, we calculate the spectrum of the two-time 
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correlation function for the induced electromotive force 
using the quantum regression theorem^ and perturba- 
tion theory. In the weak driving limit, we study how the 
Rabi splitting depends on both the strength of the driv- 
ing current which passes through the NAMR as well as 
the detuning between the frequencies of the charge qubit 
and the nanomechanical resonator. Using numerical cal- 
culations, we demonstrate that our analytical results are 
valid in the weak driving limit. Finally, we summarize 
our conclusions. 



II. MODEL 





FIG. 1: (Color online) (a) Schematic diagram of a charge 
qubit (grey loop on the left) capacitively coupled to a metal- 
lic nanomechanical resonator (NAMR) shown on the right. 
An external ac current, I{t) = lo sin{uipt) , shown in yellow, 
passes through the NAMR. Also a magnetic field B (red line 
with the up arrow) is applied to the NAMR. The magnetic 
flux $ through the SQUID is denoted by the blue upward- 
pointing arrows. The induced electromotive force V (defined 
in Eq. (Il|) at both ends of the NAMR can be used to detect 
information on the quantum-coherent coupling. A simple cir- 
cuit diagram for (a) is given in (b) . Here, Cj and Cg represent 
the capacitances for the Josephson junctions and the gate ca- 
pacitor, respectively. Also, Vg is the gate voltage applied to 
the qubit via the gate capacitor. 

As shown in Fig. [U we study a metallic NAMR, which 
is capacitively coupled to a SQUID-based Cooper pair 
box (qubit). The distributed capacitance between the 
NAMR and the superconducting island of the qubit is 
denoted by Cn{x). Here no voltage is applied to this 



distributed capacitor, in contrast to previous work^. The 
electromagnetic force drives the NAMR to oscillate and 
the induced electromotive force between the two ends of 
the NAMR can be described by 



dt M 



(1) 



Here I and M denote the length and the mass of the 
NAMR, respectively. Also, {dx/dt) = (p/M) denotes 
the velocity of the NAMR; p represents the momentum 
for the center of mass of the NAMR; x denotes the small 
displacement of the NAMR around the coordinate-axis 
origin, and x — when no ac current passes through 
the NAMR. The Hamiltonian i?NAMR for the ac-current- 
driven NAMR is given by 



Hnamr 



P 

2M 



1 



Ain^x^ - lBI{t)x, 



(2) 



where the canonical coordinate x and momentum p for 
the NAMR are assumed to satisfy the commutation rela- 
tion [x,p] = i. Hereafter, we assume h — 1. The param- 
eter n denotes the oscillating frequency of the NAMR. 
The X and p of the NAMR can be represented by the 
annihilation a and creation operators as, 
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(3) 



(4) 



Thus the Hamiltonian in Eq. ^ can be rewritten as 



IB 



V2Am 



lit). 



(5) 



Let us now assume that the SQUID has two identi- 
cal Josephson junctions, each with the Josephson energy 
Ejo and capacitance Cj. A control gate voltage Vg is ap- 
plied to the Cooper-pair box via the gate capacitor with 
the capacitance Cg. The Hamiltonian of the box can be 
writteiii^i^ai^ as 



2e^ 



{n — Ug) +Ej cosip 



(6) 



with the controllable effective Josephson energy Ej — 
2Ejo cos(7r$/$o) of the SQUID. Here, $ is the magnetic 
flux through the SQUID loop and $0 is the flux quantum. 
The total capacitance C^{x) connected to the supercon- 
ducting island is given by 

Cj:ix) = 2Cj + Cg + Cnix). 

The effective Cooper pair number Ug in the supercon- 
ducting island is Ug — CgVg/{2e). We assume that the 
distance d between the NAMR and the superconduct- 
ing island is much larger than the amplitude x of the 
oscillation of the NAMR, i.e., d ^ x. In this case, the 
distributed capacitance can be approximately written as 



C„(x)~C„(l-^) 



(7) 
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to first order in x/d. Tiie island's charging energy is 
Eq = e^/ [2Ci;(0)] and Ug is near the degeneracy point 
1/2. In these conditions, the box can be reduced to a two- 
level quantum system and the Hamiltonian in Eq. ^ can 
be reduced to 



Hr. = AEn n 



1 



C l^g 



1 



1 



in spin-1/2 notation with the quasi-spin operators 

= |0>,,(0|-|1>,,(1|, 
= |0),,(1| + |1>,,(0|, 

which are defined in the basis of the charge states |0)^ 
and Equation ([5]) is used to describe the interaction 
between the NAMR and the charge qubit. 

We now reconstruct a set of spin operators {ctz, a±} 
with 



a+|0) = |l), a_|l) = 



and 



|1) = cos-|0),-sin-|l)^, 

|0) = sin^|0), + cos^|l),, 

where the mixing angle 9 is determined by 

Ej 



(9) 

(10) 
(11) 



tan 9 



AEc {2ng - 1) ■ 



the total Hamiltonian in Eq. ((T4)) is converted into an 
effective time-independent Hamiltonian 

_ffoff = AaCr+(T_ + g (a(T+ + a) aS^ + Aa^a — ^ (a + a^) . 

(16) 

Here, the fast oscillating terms eyi'p{iujdt) and ex.'p{2iu!dt) 
have been neglected, and we also used the identity, Uz = 
2(T-|_(7_ — /, where / is the unit operator. In Eq. (|16p . the 
detuning Aa between the frequencies of the qubit and 
the ac driving current is 



Aq — Aqubit — current 



(17) 



The detuning A between the frequencies of the NAMR 
and the ac driving current is 



A EE A 



NAMR-curront 



(18) 



The interaction strength g (between the qubit and the 
NAMR) is 



9 



ffqubit-NAMR 



_ 1\ 1 C„ . 



which can be switched off when Ug — 1/2. The coupling 
strength (between the NAMR and the ac driving current) 
is 



C = ^NAMR-currcnt 



IBIo 



1 



2 V2Am 



(19) 



Note that the coupling g is proportional to the gate volt- 
age, while the other coupling strength ^ is proportional 
to BIq. 



In the new basis |0) and |1), the Hamiltonian in Eq. 
becomes 



iEc Cn 



2 dCs(O) 
with the qubit frequency 



1 



(cos 9az 



IQEl {2ng 



sin 6*172: ) X 
(12) 

(13) 



Notice that the coupling between the coordinate x of the 
NAMR and the qubit has a term proportional to the gate 
voltage (oc Ug). Thus, the gate voltage Vg can control this 
coupling. Using Eq. (fT2|) and also considering the driven 
NAMR, we can now write down the total Hamiltonian of 
the driven NAMR interacting with the charge qubit 



H = H, + H, 



NAMR- 



(14) 



Here, the Hamiltonians i?NAMR and He are given by 
Eqs. ^ and (fT^ . respectively. 

In the rotating reference frame at the driven frequency 
ujp, for both the qubit and the NAMR, through the uni- 
tary transformation 



U — exp 



(15) 



III. CORRELATION SPECTRUM OF THE 
INDUCED ELECTROMOTIVE FORCE 

The first three terms of the right hand side of iJcff 
in Eq. (jl6p describe the Jaynes-Cummings Hamiltonian, 
which was extensively studied in cavity QED. This QED 
analogue of the qubit-NAMR, described in Eq. p6|) , can 
be studied via the correlation spectrum S'y(tj) of the in- 
duced electromotive force 



V = iBl\ (a^ - a), 



(20) 



which is obtained from Eq. ^ by replacing the momen- 
tum operator p with Eq. ([?]). The correlation spectrum 
Syiijj) of the induced electromotive force V can be cal- 
culated via 

Sv{uj) = -Ref dre''^^ (F(O)V^(r)) . (21) 



Equation ([^0]) shows that the two-time correlation func- 
tion {V{0)V{t)) in Eq. ([2T|) can be calculated as 

{V{0)V{t)) oc (a(0)at(r)) + (at(0)a(T)) 

-(a(0)a(r))-(at(0)at(r)). (22) 
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A. Master equation and solutions 

To obtain the correlation spectrum, we start from the 
master equation^l of the reduced density matrix p for the 
qubit-NAMR system 



—i [Hcs^ p\ + n {2apa) — a)ap — pa'aj + 

- (2cr_p(T+ - cr+cr-p - pcr+(T_) , 



(23) 



where the latter two terms describe the decays of the 
NAMR and the charge qubit, respectively. The param- 
eters K and 7 denote the decay rates of the NAMR and 
the qubit, respectively. We also use the Markov approxi- 
mation when Eq. (|23|) is derived. For convenience below, 
we now define the number operator, N = a^cr^ -\- a^a, 
to characterize the total excitation of the qubit-NAMR. 
Obviously, N satisfies 

N\3,k) = (j + fc) |j,fc). 

Here, the index j represents the states of the charge 
qubit. When the qubit is in an excited state, we take 
j = 1, otherwise j = 0. Also k denotes the phonon num- 
ber of the oscillating NAMR, i.e., a^a\k) = k\k). 

We are only interested in the weak driving limit, i.e.. 



= ^NAMR-currc 



0. 



In this limit, we only need to consider the zero- and one- 
particle excitations; then N satisfies the condition 



N = j + k = Q, 1. 



(24) 



The Hilbert space for the reduced density matrix is now 
limited to a smaller subspace with a truncated basis 



{\3,k), j + k = 0, 1}, 



(25) 



Therefore, in this truncated basis, the density matrix el- 
ements satisfy the following equations 



dpi 



00,00 



dr 

dpoa.oi 
dr 

dpQQ,10 

dr 
dpoifii 
dr 

dpoi,io 
dr 

dpiosQ 
dr 



i^Poim - *^Poo,oi + 2kpoi,oi + 7Pio,io, 

(zA - k) Poo, 01 + igpoo,io 

+iC (Poi,oi - Poo.oo) , 

(^iAa - Poo, 10 + «5Poo,oi + i^Poi.io, 

-2kpoi,oi + ig (Poi.io - Pio,oi) 

(poo.oi — Poi,oo) , 
(^iAa - iA - K - poi,io 
+«3 (Poi,oi - Pio,io) + «Cpoo,io, 
-7Pio,io + ig (pio,oi - Poi,io) • (26) 



The other non-diagonal matrix elements Poi,oo(''')i 
Pio,oo(''")j a-iid pio,oi(''") can be easily obtained by tak- 
ing the complex conjugates of poo,oi (■''), Poo,io('''), and 



Poi,io(t), e.g., Poi,oo(t) = [poo,oi(t)]* when their solu- 
tions, e.g., Poo,oi(''")) are obtained using Eq. ^2E^ . 

In the weak driving limit, ^ <C 5, we take the pop- 
ulation in the ground state as poo,oo — 1- We also find 
that two diagonal matrix elements (pio,io and poi.oi) and 
two off-diagonal matrix elements (poi,io and pio.oi) are 
proportional to S,^. The other ones are proportional to 
^. Using perturbation theory, we only keep the terms to 
first order in ^ for the reduced matrix elements in Eq. (|26|) 
and then we can obtain 

Poo,oi = (iA - k) poo,oi + «5'Poo,io - iC, (27) 
Poo, 10 = («Aa - poo,io + «.gpoo,oi- (28) 



Applying the Laplace transformation to Eq. ([27]) and 
Eq. the solutions of the matrix element poo,oi(''') 

can be easily obtained as 



Poo,oi(t) = 77126 
with parameters 



(29) 



A1A2 

PmnPOO,Ol(0) -I- XmnPOO,lo(0) + i^I^n 



Other parameters Pmn, Xmn, '^mn, Ai, and A2 in Eq. ([29 
are 



Xr, 



Pmn 



{iAa 



(Am ~ A„ 

A,„ + (iAa - 



(Arn An) 

g 



i (Am — A„) 



(30) 



and 

Am. = 



{-l)"^^{5~^(n~l)) +452-(A,-fA) 



for m{^n) = 1, 2. Where we define the frequency de- 
tuning 



(5 = Aa - A = tj„ - f], 
and the parameter F is given by 

F = 



7 
4' 



(31) 
(32) 



The parameters can be further expressed as, A^^ 
-I- iifm , with real part 



F. 



r + ^(-i)^ 



2\ 4 





'1 , 




sin 


— arctan 






2 \ 





and imaginary part 
fni = ^ |(-1)™ {a^ 



2_ 


'1 / 




* COS 


— arctan 






2 \ 





- [Aa + A) 
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Here, the parameters a and h are 



with 



b = 26 [k- - 
2 



B. Correlation spectrum 

The correlation function, e.g., (a^^(O)a(r)), is given by 
(at(0)a(r))=Tr{a(OM(r)} (33) 

with 

A{t) = U{T)p{id)a\0)U\T). (34) 

Using the quantum regression theorem^l'^ the correlation 
function in Eq. (|33p can be written as 



(at(O)a(T)) -^oi,oo(r), 



(35) 



with 



Aoi,oo(t) = i?i,i2e-^^" + i?i,2ie"^^" + £*■ (36) 
Here, the initial operator A{0) is assumed to be 

A(0) =p""at(o) (37) 

when we calculate the time-dependent matrix element 
^01,00 (■'') in Eq. ([36]) . The "ss" in the superscript of the 
reduced density matrix p denotes the "steady state" . The 
parameters Ha,i2 and Ha,21 in Eq- P6p are expressed as 

HA.mn = A^m„^01,0o(0) + Xmn^lO,Oo(0) " «Ci^m„, (38) 

with the subscript either mn = 12 or mn = 21. The pa- 
rameters Aoo.oi(O) and Aoo.io(O) denote the matrix ele- 
ments of the operator A{{)) in the truncated basis defined 
in Eq. e.g., 

^oo,oi(0) = (00|A(0)|01) = (00|p^^at(o)|oi). (39) 

These matrix elements ^oi,oo(0) and Aio,oo(0) can be 
straightforwardly obtained as 



^oi,oo(0) - Poi,01> 
^io,oo(0) = PiQfiii 



(40) 
(41) 



where and pfg denote the "steady-state" matrix 

elements of the reduced density matrix p. 

Using the same procedure, other correlation functions 
can also be obtained as 

(a(0)a+(T)) = Tr{a+(0)i3(r)} =i3oo,oi(T), (42) 
(a(O)a(r)) = Tr {a(0)C(r)} = Coi,oo(t), (43) 
{a\Q)a\T)) = Tr{a+(0)7^(r)} =7^oo,oi(r), (44) 



B{t) = U{T)p{0)a{0)U\r), 
C{t) = [/(r)p(0)a(0)t/t(r), 
D{t) = [/(T)p(0)at(0)t/t(r), 



(45) 
(46) 
(47) 



and 



Soo,oi(t) = HB,i2e-^'^ +HBaie-^'^ +e, (48) 
C^oi,oo(t) = il£,i2e-^i^+i?2/,2ie"^^'"+e*, (49) 
DooM^) = HDA2e-^'^ + HD^2ie-^'^ +e. (50) 
Here, the parameters HB.mn, Hc,mm and Hu^mn are 

Hc\,nn = A^m„C'oi,Oo(0) + Xm«C'lO,Oo(0) - iS.'^^n^ (51) 

and 

Hx.rnn — P'mnXQQ^Qi{Q) + Xmn-'^OO.IO (0) -|- i^Vrnn- (52) 

Here, the subscript X can be either B or D. With 
the same meaning as in Eq. (j39p . the parameters, e.g., 
C'oi,oo(0)j represent the matrix elements of the operators 
C(d), B(0), and 1^(0) in the truncated basis in Eq. ((25|l . 
We note that the initial conditions 

B{Q) = C(0) =p"'*a(0), (53) 
D{Q) = A(0) =p^"a^(0) (54) 

are used when Eqs. (jHHHI) are derived. The matrix ele- 
ments, e.g., i?oo.oi(0), can also be obtained as 

Soo,oi(0) = pgg^oo. (55) 

using the quantum regression theorem with the matrix 
element Pqq qq of the reduced density matrix p. 

By using Eq. and Eqs. (|42p4)) . the two-time corre- 
lation function in Eq. ([22|) for the induced electromotive 
force is given by 

(F(O)U(t)) (X Aoi^oo(t) + Soo,oi(t) 

-C'oi,oo(t) - Doo,oi(t)- (56) 

Based on the above results, the two-time correlation func- 
tion in Eq. ()56p is further simplified to 



(F(O)U(t)) cx m2e"^^"+Ai2ie-^^" 



+/ie"^i" + he 



-\% 7 



(57) 



with 



/i 

/2 



A'l2Poi,01 + Xl2PlO,01) 

A^2iPoloi + X2lPlO,01- 

Then, replacing (U(O)U(t)) in Eq. ^ by Eq. ([57| . 
and integrating, the spectrum S'y(w) in Eq. (|2ip can be 
expressed as 

Sy{u) « 444^ Re 



2Af 

B2/2f7riRe(^i2) 



Ai2ie 



ipi)lu\(pi 



2M {u;~^iY+Tl 
BH^n r2Re(^2i) - (w - V2) Im(Ai2i) 



2M 



(w - (^2)' 



(58) 
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In Eq. ([58)) . we have neglected the terms proportional 
to the amplitudes /i and f2- This because the ratios, 
e.g., {fi/nu), are proportional to which is neghgibly 
small in the weak driving limit. In this case, we need only 
consider the two leading terms, which are proportional to 
/ii2 and /i2i, as shown in Eq. ([58]) . 

As shown in Figs.[2l[3l[4l there are two dominant peaks 
in the S'v'(w) spectrum. The distance (splitting frequency 
Auj) between these two peaks is 



"1 




— arctan 




2 





(59) 



which is determined by the frequency detuning d and 
the decay rates k and 7 for the NAMR and the charge 
qubit. Then, the information of the coherent coupling 
between the charge qubit and the NAMR can be obtained 
by Eq. HMD. 
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FIG. 3: The spectrum Sy(a;) versus ui for three different val- 
ues of the detuning S, e.g., 5 — 0, 0.4, 0.8. The detuning 5 is 
defined in Eq. (|3ip . The total excitation number A'^ is taken 
here as A = 1. Other parameters and units are the same as 
in Fig. H 



CO 



50 


50 


50 




^=0.02 



^= 0.03 



1 



^= 0.04 



1 



-1.5 



-1.0 
CD 



-0.5 



FIG. 2: The spectrum Sv(i^) versus lu for three different val- 
ues of the total number of excitations A, e.g., A = I, 2, 3. 
We take the values for other parameters as 5 — 0.2, g — 0.2, 
^ = 0.02, K = 0.004 and 7 = 0.004. Here, aU these parameters 
are in units of I GHz. 



FIG. 4: The spectrum S'v(i^) versus lu for three different val- 
ues of the driving strength 5, e.g., £, = 0.02, 0.03, 0.04. The 
coupling strength ^ = ^NAMR-currcnt is defined in Eq. (|19|l . 
The total excitation number N is taken here as A = 1. Other 
parameters and units are the same as in Fig. [2] 



C. Numerical results 

To test the validity of our analytical calculations ob- 
tained using both quantum regression theorem and per- 
turbation theory, we now study the correlation spectrum 
Svi(^) numerically. 

In Fig. [21 the spectrum Sv{to) versus frequency w is 
plotted with parameters S = 0.2, g — 0.2, ^ — 0.02, 
K = 0.004, and 7 — 0.004, for different total excitation 
numbers A^, e.g., A^ = 1, 2, 3. Here, we take 1 GHz as 
the unit for all these parameters. 

Figure[2]shows that: (i) there are two prominent peaks, 
which means that the approximation in Eq. (|58p is valid 
in the weak driving limit; (ii) the increase of the total 
number A^ of the excitation does not obviously change the 
heights and the splitting frequency Au of the two leading 
peaks. Therefore, Fig. [2] verifies that the approximation 



with the truncated basis (A^ = 1) is valid in the limit of 
weak driving <C 5). 

Now in the truncated basis (A^ = 1), let us demonstrate 
the effects of the frequency detuning S [see Eq. (PTjl ] and 
the driving strength ^ [see Eq. (|19p ] on the splitting fre- 
quency Auj. In Fig. [3l the spectrum S'v(tf) versus fre- 
quency u is plotted for different detunings 6, and the 
same other parameters as in Fig. [2j Figure [3] shows that 
the detuning S affects both the heights and the frequency 
splitting Au! of the two peaks. A larger S corresponds to 
a larger frequency splitting Alu. When S is increased, 
the heights of the two peaks change asymmetrically, i.e., 
one peak becomes higher than another one when increas- 
ing the detuning S. (Similar results were found in Ref.^). 
This means that our probing approach does not work well 
when the frequency detuning (S = 0.8) is much larger 
than the interaction strength {g = 0.2). In practice, we 
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should assume that the frequency detuning S and the in- 
teraction strength g is of the same order. Also, Fig. [3] 
shows that the detuning i5 does not affect the number of 
peaks, which means that the approximation in Eq. (j58p 
is valid in the weak driving limit. 

Similarly, in Fig. HI the spectrum Sv(j-^) versus fre- 
quency u is plotted for different driving strengths ^ and 
the same other parameters as in Fig. [21 when the total 
excitation number = 1. Figure [4] shows that a weak 
driving strength ^ does not significantly affect the fre- 
quency splitting Auj of the two peaks, and also it does 
not affect the number of peaks. Moreover, changing the 
driving strength ^ does not affect the heights of both 
peaks. 

IV. CONCLUSIONS 

We have proposed an "intrinsic probing" approach 
to demonstrate the coherent coupling between a driven 
metallic NAMR and a charge qubit. This metallic 
NAMR can act as an intrinsic detector when a weak driv- 
ing current passes through it. Using the quantum regres- 
sion theorem and perturbation theory, we have calculated 
the correlation spectrum of the electromotive force be- 
tween two ends of the NAMR. This spectrum can be used 
to demonstrate QED analogues in the NAMR-qubit sys- 



tem, e.g., the vacuum Rabi splitting related to the coher- 
ent coupling strength of the charge qubit to the NAMR. 
The numerical calculations confirm the validity of the 
analytical results. 

In our proposal, no additional measurement instru- 
ments needs to be integrated, in contrast to the proposal 
in^. The NAMR acts as both the cavity and the detec- 
tor. Therefore, it is easier to be fabricated. Our proposal 
can also be generalized to the case where many qubits 
are coupled to a NAMR. In this case, the information of 
many qubits can also be readout via the spectrum of the 
electromotive force. Recent experiments^'^ indicate that 
our proposal is experimentally realizable. 
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